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Abstract 

We have studied the e + e~ — ► coir cross section in the y/s interval 1000-1030 MeV 
using the 7r + 7r~7r 7r and 7r 7r°7 final states with a sample of ~ 600 pb _1 collected 
with the KLOE detector at DA<J>NE. By fitting the observed interference pattern 
around M^ for both final states, we extract the ratio of the decay widths T(lo — > 
7r°7)/r(w — > 7T + 7T~7r ) = 0.0897±0.0016 and derive the branching fractions BR(w — > 
tt+tt-tt ) = (90.24 ± 0.19)%, BR(w -» tt°j) = (8.09 ± 0.14)%. The parameters 
describing the e + e~ — ► 7r + 7r~7r°7r° reaction around M^ are also used to extract the 
branching fraction for the OZI and G-parity violating <fi — > ton decay: BR(</> — ► 
(jvr°) = (4.4 ±0.6) x 10" 5 . 

Key words: e + e~ collisions, rare 4> decays, VMD, OZI violation, Isospin violation 



1 Introduction 



At low energy, below 1.4 GeV, the e + e~ — > uir° cross section is largely dom- 
inated by the non-resonant processes e + e _ — ► p/p' — > ujtt . However, in the 
region around M^, a contribution from the OZI and G-parity violating de- 
cay (f)^uJ7Y is expected. This strongly suppressed decay (O(10 -5 )) can be 
observed via interference with the non-resonant process, showing up as a dip 
in the total cross section dependence from y/s. 
The e + e~ — ► ujtt cross section as a function of ^/s is parametrized as [1]: 



Cr(y/s) = CTnr(Vs) 



„ M^T„ 

Zj - 



n*W*) 



where a nr {^/s) is the bare cross section for the non-resonant process, Z is 
the interference parameter (i.e. the ratio between the decay and the non- 
resonant process amplitudes), while M^, T^ and D<p are the mass, the width 
and the inverse propagator of the meson respectively. The non-resonant 
cross section, in the energy range of interest, increases linearly with y^i. A 
model independent parametrization is used to describe the non resonant part: 

In this work we study two different final states: 7r + 7r _ 7r°7r° and 7r°7r°7, corre- 
sponding to the uj decay in 7r + 7r~7r° ad 7r°7, respectively. From the 7r + 7r~7r 7r° 
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analysis we have measured the Z parameter that has been used to determine 
the 0— >um° branching ratio (BR). In the case of 7r°7r°7 we expect contributions 
also from — > pir and — > S7 intermediate states, being S a scalar meson. 
In [2] we have shown that at y/s ~ M^ the interference between — »• £7 and 
e + e~ — > urn events, evaluated by fitting the M^^-M^ Dalitz plot, is small. 
Therefore a fit to the cross section interference pattern for 7r°7r°7 final state 
will provide information about the e + e" — > p/p' —>■ urn process and the reso- 
nant decays (p^un and — > /wr . Moreover, combining the results we extract 
the ratio r(u; — > 7r°7)/r(u; — > 7i + ir~7i°). 



2 The KLOE detector 



The KLOE experiment operates at DA$NE [3], the Frascati 0-factory. DA$NE 
is an e + e~ collider running at a center of mass energy of ~ 1020 MeV, the 
mass of the meson. Equal energy positron and electron beams collide at an 
angle of 7T-25 mrad, producing mesons nearly at rest. 

The KLOE detector consists of a large cylindrical drift chamber, DC, sur- 
rounded by a lead-scintillating fiber electromagnetic calorimeter, EMC. A su- 
perconducting coil around the EMC provides a 0.52 T field. The drift cham- 
ber [1], 4 m in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense 
wires and 37,746 aluminum field wires. The chamber shell is made of carbon 
fiber-epoxy composite and the gas used is a 90% helium, 10% isobutane mix- 
ture. These features maximize transparency to photons and reduce Kl —>■ Ks 
regeneration and multiple scattering. The position resolutions are a^^lSO pm 
and <r 2 ~ 2 mm. The momentum resolution is a(p±)/p± ~ 0.4%. Vertices are re- 
constructed with a spatial resolution of ~3 mm. The calorimeter [5] is divided 
into a barrel and two endcaps, for a total of 88 modules, and covers 98% of the 
solid angle. The modules are read out at both ends by photo-multipliers, both 
in amplitude and time. The readout granularity is ~ (4.4x4.4) cm 2 , for a total 
of 2440 cells arranged in five layers. The energy deposits are obtained from 
the signal amplitude while the arrival times and the particles positions are ob- 
tained from the time differences. Cells close in time and space are grouped into 
calorimeter clusters. The cluster energy E is the sum of the cell energies. The 
cluster time T and position R are energy-weighed averages. Energy and time 



resolutions are a E /E = 5.7%/ JE (GeV) and a t = 57 ps/JE (GeV) © 100 ps, 
respectively. The KLOE trigger [B] uses both calorimeter and chamber infor- 
mation. In this analysis the events are selected by the calorimeter trigger, 
requiring two energy deposits with E > 50 MeV for the barrel and E > 150 
MeV for the endcaps. A cosmic veto rejects events with at least two energy 
deposits above 30 MeV in the outermost calorimeter layer. Data are then an- 
alyzed by an event classification filter [7], which streams various categories of 



events in different output files. 



3 Data analysis 



All available statistics collected at the <p peak during 2001-2002 data-taking 
periods, corresponding to 450 pb" 1 , have been analyzed. Moreover four scan 
points (at 1010, 1018, 1023 and 1030 MeV) of ~10 pb" 1 each and the off-peak 
run at 1000 MeV of ~100 pb _1 acquired in 2005-2006 have been included 
in this analysis. The luminosity is measured with 0.5% absolute precision 
counting large angle Bhabha scattering events [8]. Data taken at the peak 
are grouped in center of mass energy bins of 100 keV width. For all the other 
points, close y/s values are grouped together and the average center of mass 
energy is evaluated by weighting with luminosity. 

In the following, the visible cross section for a given final state (j =4ir ,7nr , -f) 
is defined a J vis = Nj/ejL[ nt , where Nj, Ej and L int , are the number of signal 
events, the analysis efficiency and the luminosity, respectively. 
The visible cross section is related to the bare cross section through the radi- 
ator function H (see Sec. @J: 

<4(v^ = Jds'H(s, s')^(s') = ffi(Vi) * SUsft- 



3.1 e + e~ — > utv° —> 7r + 7r~7r°7r° 



In the 7r + 7r~7r°7T analysis, data are filtered by selecting events with the ex- 
pected final state signature: two tracks with opposite curvature connected to 
a vertex inside a small cylindrical fiducial volume (p < 4 cm and \z\ < 6 
cm) around the Interaction Point (IP) and four neutral clusters in the prompt 
Time Window (TW), defined as |T 7 — i? 7 /c| < MIN(4o"(,2 ns). To minimize 
contamination from machine background, we require for the clusters a min- 
imal energy of 10 MeV and an angle with respect to the beam line in the 
interval 22° < # 7 < 158°. 

With this selection, at <fi resonance, the main background contributions come 
from <fi — > K$Kl — > 7r + 7r _ 7r 7r /7r 7r 7r/iz/ and <fi —>■ K + K~ , with K ± — > ^n , 
which have the same final state signature. Other two resonant background 
components (0 — > 777 with 77 — > ir + TT~7r°, and — > 7T + 7r~7r°) mimic the final 
state signature because of additional clusters due to accidental coincidence of 
machine background events and/or shower fragments (cluster splitting). An 
additional non-resonant background contribution of the order of few percent 



from e + e — > tt + tt 7i°ir , dominated by the ai(1260)7r intermediate state, is 
also expected for all y/s values [9] . 

A global kinematic fit (iXiof = 8), imposing total four-momentum conservation 
and proper time of flight (TOF) for photons coming from the charged vertex, 
improves both the signal/background separation and the determination of 
the photon energies. The resulting \ 2 (Xxfit) * s use d to select signal enriched 
(Xxfit < 50), S cvts , and background dominated (xxfit > 50), B evts , samples. In 
S cvts sample the overall contamination from resonant background at y/s ~ M^ 
is about 12% while becomes negligible outside <fi resonance. The contribution 
of the ai(1260)7r background is about 4% for all y/s values. The signal selection 
efficiency in the S ev ts sample is evaluated by Monte Carlo (MC) and corrected 
with data/MC ratios for tracking, vertexing and clustering. The resulting value 
£4tt ~ 38% is dominated by the selection requirements and shows a small 
dependence on y/s, which is taken into account in the evaluation of a^/ 



VIS' 



The signal counting is performed for each y/s bin by fitting the tt° recoil 
mass (M rcc ) distribution for both, S evts and B cvts samples with MC signal and 
resonant, non-resonant background shapes. The fit procedure is based on a 
likelihood function which takes into account both data and MC statistics. 
In Fig. [TJa-d, data-MC comparisons for events in the most populated bin of 
y/s are shown. In Fig. dJe-f, the M rec distribution is also shown for the two 
outermost center of mass energies. 

The results are summarized in Tab.[IJ where the signal counting (N An ) and the 
visible cross section (cr^ s ) are reported for all y/s bins. Errors on a^ s include 
statistics, background subtraction and a systematic error of 0.75%, dominated 
by the luminosity measurement (0.5%). The other systematics being due to 
cosmic ray veto, event counting and final state radiation. 

We have checked the stability of the result by repeating the whole analysis 
chain with a wide variation of the selection criteria on: (a) minimum cluster 
energy and angle in the TW definition; (b) value of Xxat cu ^ used to select S cvts 
and B evts samples; (c) distribution used for the B evts sample when performing 
signal counting; (d) data/MC efficiency correction curves. The resulting a^f s 
values are then used to estimate the systematic error of the measurement. 
The fit to the visible cross section as a function of y/s described in Sec. H] 
is repeated for all of these variations. The quoted systematics uncertainty is 
calculated as the quadratic sum of RMS's obtained from variations (a)-(d). 

3.2 e + e~ -+ un° -> 7r°7r°7 



The selection for 7r°7r°7 starts requiring five neutral clusters in the prompt 
Time Window with E 1 > 7 MeV and polar angle | cos6* 7 | < 0.92. After a first 
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Fig. 1. Data-MC comparison for 7r + 7r~7r°7r° events taken at \fs = 1019.75 MeV: 
X 2 of kinematic fit for S cvts (a) and B cvts (b) samples, angle between charged pions 
in the uj rest frame (c) and 7r° recoil mass (d) for S cv t s events. The last distribution 
is shown also at y/s = 1000.10 MeV (e) and y/s = 1029.95 MeV (f). 

kinematic fit (Fitl, iVdof = 9) imposing total 4-momentum conservation and 
time of flight, photons are paired to 7r°'s by minimizing a \ 2 function, built 
using the invariant mass of the two 77 pairs. A second kinematic fit (Fit2, 
iVd f = 11) imposes also the constraint on n° masses. 



The background is rejected by requiring Xp it2 /Ndf < 5 and AM, 



77 



|M 77 - 



M n \ < 5 a 77 , where M 77 and cr 77 are evaluated using the photon momenta from 
Fitl. After these cuts the remaining sample is dominated by e + e~ —>■ uir — > 
7r°7r°7 and (ft — > S7 — > 7r 7r°7 events. Signal is then selected neglecting the in- 
terference between these two processes and cutting on the intermediate state 
mass. Defining M 7r7 as the closest mass to M^ of the two 7r°7 combinations, 
only events satisfying the requirement 750 < M ny < 830 MeV are retained. 
The residual background contamination (~ 20% at the <p peak) comes pre- 
dominantly from — > 777 — > 7r 7r°7r°7 events, where two photons are lost or 
merged. 
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Fig. 2. Data-MC comparison for 7r°7r°7 events taken at y/s = 1019.75 MeV: (a) 
normalized x 2 of the second kinematic fit after acceptance cuts (Fit2), (b) normal- 
ized two-photon invariant mass after Xp it2 cut, (c) 7r°7 invariant mass after AALy 
cut, (d) costp distribution after all analysis cuts; M n7 distribution at y/s = 1000.10 
MeV (e) and y/s = 1029.95 MeV (f). 

In Fig. [2ja-d, data-MC comparison for events in the most populated y/s bin 
are shown. The if) variable is the minimum angle between the photon and 
the 7r°'s in the di-pion rest frame. A good agreement is observed both after 
acceptance selection and after applying analysis cuts. The comparison for the 
M 7f7 distribution is also shown for the two outermost center of mass energies, 
where the background contribution is practically negligible (Fig. [2je-f). 



The overall signal selection efficiency is evaluated by applying the whole anal- 
ysis chain to signal MC events: £7^7 ~ 46%, almost independent from y/s. 
The value obtained for each bin, together with the corresponding integrated 
luminosity, has been applied to the signal counting to obtain the visible cross 
section (<r vis ' ). Results are summarized in Tab. [TJ Errors on cr vis include 
statistics, background subtraction and a systematic error of 0.6%, dominated 
by the luminosity measurement (0.5%). The other systematic effects, related 



to cosmic ray veto and to the event classification filter, have been evaluated 
with downscaled data samples. 

The stability of the results has been checked by repeating the fit with: (a) a 
variation of the selection criteria on the TW definition and on xlit2 an d M ny 
cuts; (b) a rescaling of the background components according to the value 
obtained by fitting background-enriched distributions on data with the cor- 
responding MC components; (c) a subtraction to the visible cross section of 
the interference effect between <p —> S'-f and the signal. This interference con- 
tribution is extracted from a fit to the Dalitz plot of the 7r°7r°7 final state [2] 
at the 4> peak, which provides the parameters describing the two processes. 
These values are then used to evaluate for each yj~s the interference contribu- 
tion, which is subtracted to the corresponding a^ 7 . The largest contribution 
is 1.4%. 
Table 1 

Signal counting, visible cross section and radiative correction for e + e~ — > 7r + 7r - 7r°7r° 
and e + e~ — > 7r°7r°7 events. The radiative correction values are calculated using for 
the bare cross section parameters in Tab. [2j 



v^ (MeV) 


N 4,r 




a An (nb) 

VIS V ' 


rati, 


N^t 


* n V (nb) 

VIS \ ' 


raa 


1000.10 


221917 ± 


562 


5.72 ±0.05 


0.885 


27110 ± 167 


0.550 ±0.005 


0.889 


1009.90 


25968 ± 


173 


6.20 ±0.06 


0.893 


2958 ± 56 


0.581 ±0.012 


0.895 


1017.20 


16209 ± 


171 


5.71 ± 0.08 


0.924 


2108 ± 46 


0.564 ±0.018 


0.923 


1018.15 


22167 ± 


158 


5.60 ±0.06 


0.933 


2557 ± 60 


0.541 ± 0.014 


0.944 


1019.30 


4799 ± 


90 


5.88 ±0.12 


0.918 


322 ± 22 


0.480 ± 0.034 


0.964 


1019.45 


58077 ± 


340 


5.89 ±0.06 


0.912 


6058 ± 101 


0.497 ±0.009 


0.962 


1019.55 


93596 ± 


445 


5.93 ±0.05 


0.908 


9516 ± 130 


0.487 ±0.008 


0.960 


1019.65 


171571 ± 


888 


5.98 ±0.05 


0.904 


18349 ± 189 


0.509 ± 0.007 


0.958 


1019.75 


326774 ± 


872 


6.04 ±0.05 


0.900 


34049 ± 282 


0.505 ±0.006 


0.955 


1019.85 


256008 ± 1248 


6.08 ±0.05 


0.896 


26124 ± 234 


0.508 ± 0.006 


0.951 


1019.95 


35850 ± 


263 


6.20 ±0.07 


0.893 


3510 ± 74 


0.491 ±0.011 


0.948 


1020.05 


17971 ± 


167 


6.21 ±0.08 


0.889 


1843 ± 52 


0.516 ±0.016 


0.944 


1020.15 


8190 ± 


132 


6.23 ±0.11 


0.886 


702 ± 32 


0.501 ± 0.024 


0.940 


1020.45 


9657 ± 


117 


6.41 ± 0.09 


0.878 


667 ± 31 


0.488 ± 0.024 


0.928 


1022.30 


16931 ± 


111 


7.24 ± 0.08 


0.869 


1891 ± 43 


0.612 ±0.018 


0.891 


1023.00 


29611 ± 


177 


7.41 ± 0.07 


0.871 


3101 ± 61 


0.619 ±0.013 


0.888 


1029.95 


33681 ± 


186 


7.84 ±0.07 


0.887 


3896 ± 65 


0.689 ±0.013 


0.892 



4 Fit results and uj branching ratios extraction 



The measured values of visible cross section, shown in Tab. [U are fitted with 
the parametrization (JTJ convoluted with the radiator function [10]. The free 
fit parameters of the bare cross section are: a J , $i(Zj), $f(Zj) and a'j, where 
j represents either the 4tt or 77777 final state. In Fig. [3] data points with the 
superimposed fit result are shown for both channels. The values of the pa- 
rameters are reported in Tab. [2j The second error quoted is the systematic 
uncertainty. It is evaluated as the quadratic sum of the RMS of the param- 
eters extracted by repeating the fit with different conditions, as described in 
the previous section. The resulting x 2 ,M&of are 11.79/13 (P(x 2 ) — 54%) and 



4.78/13 (P(x 2 ) = 98%) for 7r vr°7 and Tr+Tr-vrV channel, respectively. The 
fit has been repeated using a VMD model [1] based on p and p' intermediate 
states for the non-resonant term of the cross section. Results are in agreement 
with the linear parametrization within one standard deviation] 1 1 



Table 2 

Fit results for the e + e~ 



7r + 7r 7T 7T and for e + e 



7r°7r°7cross section. 



Parameter 


e + e — > 7r + 7r 7r°7r° 


e + e — ► 7r°7r°7 


a [nb] 


7.89 ±0.06 ±0.07 


0.724 ±0.010 ±0.003 


Re(Z) 


0.106 ±0.007 ±0.004 


0.011 ±0.015 ±0.006 


Qm(Z) 


-0.103 ±0.004 ±0.003 


-0.154 ±0.007 ±0.004 


a' [nb/MeV] 


0.064 ± 0.003 ± 0.001 


0.0053 ± 0.0005 ± 0.0002 



After removing common systematics on the luminosity (0.5%), from the two 
measurements we obtain: 



a (u -> tt°7) 



a [u> 



TX TX TX 



0\ 



0.0918 ±0.0016 



(2) 



Taking into account the phase space of the two decays [T], the ratio of the 
partial widths is: 



T(u -> 7T°7) 



V(u 



TX TX IX 



Q\ 



0.0897 ±0.0016 



(3) 



Since these two final states account for 98% of the uj total width, we use the 
T{uj —* 7r°7)/r(ct; —* 7r + 7r _ 7T°) ratio and the sum of rarer BR's [UJ to obtain 
from a fit: 



The fitting function used for the non-resonant term is: 



<(£) 



ifn", 



0"n 



E 3 



m 2 p U p {E) + AjmlU^E) 



P { (E) 



mfilpim^) + Ajrri^Up'im^ 



P{(m 4 



where U. p m (E) 



(m 2 



p(') 



E — iET p n\{E)) is the vector meson propagator while 



Pf(E) describes the energy dependence of the phase space volume of the final state. 
For the latter we assume the approximation of an infinitely narrow u) meson. The 
free parameters in this case are a^ and the real number Aj. This parametrization 
allows us to compare directly the value of <Tq obtained from this model with the 
value in Tab. [2j The width of the p meson has a dependence with energy T p (E) = 
Tp(mp)(mp/E)(p 7r (E)/pT T (mp)), while the p' width is assumed fixed. m p o\ and T p (/) 
values are taken from [TT]. For the parameter A we obtain: —0.15 ±0.04, where the 
error takes into account only the contributions from Tab. [TJ 
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Fig. 3. Cross section fit results for the e + e — > tt^tt 7r°7r° (top) and e + e — > 7r°7r°7 
(bottom) channels. Black dots are data, solid line is the resulting fit function. 



BR(iV -> TT + TT-TT ) 

BR(u -> vr°7) 



(90.24 ±0.19)% 
( 8.09 ±0.14)% 



(4) 
(5) 



with a correlation of -71%. Comparison between our evaluation and the values 
in [11] is shown in Fig. HI Our result for BR(w — > 7T 0/ y) is less than the PDG 
value by three standard deviations. It is in good agreement with the recent 
prediction in [T2] . 
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Fig. 4. Branching fraction for the two main oj decay channels. The black square is 
the KLOE fit result, while the black dot is the constrained fit result in [11]. The 
gray ellipses are the 68% C.L. regions. 



5 BR(0 — ► lutt ) evaluation 



The measured a^ and Z±„ parameters of the 7r + 7r 7r°7r° final state are related 



to the BRO 



BR( 



>UJ1T 



through the relation: 



U7T°) = 0-0(771^)1^4^1 



O 



4-7T I 17 I 



Ml 



o> BR(uj -> tt+tt^tt ) 12vrr(0 -> e+e") 



(6) 



where ao(m^) is the total cross section of the e + e~ — > u;7r process and 0$ 
is the peak value of the bare cross section for the <fi resonance. Using the 
parameters obtained from the 7r + 7r~7r 7r° analysis, the T ee measurement from 
KLOE [13] for the evaluation of a^, and our value for BR(u; — > tt + 7t~7t°) we 
extract: 



BR(0 



UJ7T 



(4.4 ±0.6) x 10" 



(7) 



The error is reduced by a factor of two with respect to the best previous 
measurement from the SND experiment [1], which is in agreement with our 
result. 
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6 Conclusions 



Using a sample of 600 pb" 1 collected at center of mass energy between 1000 
and 1030 MeV, we have measured the cross section parameters for the two 
processes e + e~^7r + 7r~7r°7r° and e + e~^7r°7r°7, obtaining the ratio T(uj — ► 
7r°7)/r(u; — » 7T + ir~7T°) with an accuracy of 1.8%. This ratio, together with 
the unitarity relation and the BR measurements on the other u decay chan- 
nels, substantially improves the accuracy on the dominant u branching frac- 
tions giving a value for BR(w — > 7r°7) which is three standard deviations lower 
than the Particle Data Group fit [11] . Moreover, the parameters describing the 
e + e~ —>■ 7r + 7r - 7r°7r° reaction around M^ are used to extract the most precise 
BR measurement of the OZI and G-parity violating decay (p^un . 
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